Pairing correlations in nuclei on the neutron-drip line 
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Paring correlations in weakly bound nuclei on the edge of neutron drip line is studied by using 
a three-body model. A density-dependent contact interaction is employed to calculate the ground 
state of halo nuclei 6 He and n Li, as well as a skin nucleus 24 0. Dipole excitations in these nuclei 
are also studied within the same model. We point out that the di-neutron type correlation plays 
a dominant role in the halo nuclei 6 He and 11 Li having the coupled spin of the two neutrons S=0, 
while the correlation similar to the BCS type is important in 24 O. Contributions of the spin 5=1 
and S=0 configurations are separately discussed in the low energy dipole excitations. 
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PACS numbers: 21.30.Fe,21.45. 



I. INTRODUCTION 
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It is now feasible to study the structure of nuclei on 
the edge of neutron drip line. Such nuclei are expected 
to have unique properties influenced by the large spatial 
distribution of weakly bound valence neutrons, for ex- 
amples, halo, skin, new magic numbers and strong soft 
dipole excitations. 

Two-neutron halo nuclei (sometimes referred to as bor- 
romean nuclei when there is no bound state between a 
valence neutron and a core nucleus) like 6 He and n Li 
have been often described as three-body systems consist- 
ing of two valence neutrons interacting with each other, 
and with the core [1-7]. The three-body Hamiltonian 
with realistic two-body interactions has been solved by 
the Faddeev method [6, 7]. On the other hand, Bertsch 
and Esbensen have developed a three-body model with a 
density-dependent delta interaction among the valence 
neutrons [1]. They have subsequently extended their 
model by taking into account the effect of the recoil 
of the core nucleus [3]. They showed that the density- 
dependent contact force well reproduces the results of 
Faddeev calculations even though the radial dependence 
of the adopted interactions are quite different [3] . To the 
date, the most sophisticated many-body calculations of 
light nuclei include also three-body forces which play an 
important role in obtaining the correct binding energies 
of light nuclei [8]. To a large extent, such three-body 
forces can also be simulated effectively by a density de- 
pendent force. 

Recently, a Hartree-Fock-Bogoliubov (HFB) model has 
been applied to study a di-neutron structure of the drip 
line nuclei. The dipole response has also been stud- 
ied with the quasi-particle random phase approxima- 
tion (QRPA) taking into account the continuum effect 
[9]. Cluster models have also been often used to study 
the borromean nuclei including the dipole excitations 
[10, 11]. 

In this paper, we undertake a detailed discussion on the 
ground state as well as the dipole response of neutron- 
rich nuclei using a three-body model with the density- 
dependent delta force, paying special attention to the di- 



neutron structure of valence neutrons. We particularly 
study the borromean nuclei, 6 He and 11 Li, and also an- 
other drip line nucleus 24 as a comparison. Since 23 
is bound, 24 is not a borromean nucleus by definition, 
although the root-mean-square (rms) radius indicates a 
feature of very extended neutron wave functions. The 
model we use is essentially the same as that of Bertsch 
and Esbensen [1-3], while the interaction is adjusted to 
fit the separation energy of each drip line nucleus. 

The paper is organized as follows. In Sec. II, we dis- 
cuss the three-body model and the adopted two-body 
interactions. The results for the borromean nuclei are 
compared with those for the drip line nucleus 24 in Sec. 
III. A summary is given in Sec. IV. 



II. THREE-BODY MODEL 

We consider a three-body system consisting of two va- 
lence neutrons, and an inert core nucleus with the mass 
number A c . We use the same three-body Hamiltonian as 
in Ref. [3]. That is, 



H = h nC (l) + h nC (2) + V nn + 



P1P2 

A r m 



(1) 



Here, h n c is the single-particle Hamiltonian for a valence 
neutron interacting with the core and is given by 



P 



hnc = ^- + V n c(r), 



(2) 



where /1 = mA c /(A c + 1) is the reduced mass. The re- 
duced mass \i, together with the last term in Eq. (1), 
originate from the recoil kinetic energy of the core [3]. V nn 
is the interaction between the valence neutrons given by 



V nn (r 1 ,r 2 ) = S(n - r 2 ) v + 



1 + exp[(n - R P )/a p ] 



(3) 

It is well known that the delta force (3) must be sup- 
plemented with an energy cutoff E cut in the two-particle 
spectrum. In terms of the energy cutoff E cut and the 



2 



scattering length a nn for nn scattering, the strength for 
the delta interaction v is given by [3] 



2v 2 h 2 



v 



m 



2a n 



2k n a n 



(4) 



where E cut — h 2 k 2 /m. The parameters for the density 
dependent part, i.e., v p , R pi and a p , are adjusted in order 
to reproduce the known ground state properties for each 
nucleus. We specify the value of the parameters below. 

We diagonalize the Hamiltonian (1) in the model space 
of the two-particle states with the energy ei + e 2 < 
(A c + l)E cu t/A c [3], where e is a single-particle energy 
of valence particle. We use a Woods-Saxon potential for 
V n c to generate the single-particle basis: 



Vnc(r) = Vb ( 1 - 0.44/ so r5(/ ■ s)- — 

r-R 



x 



1 + cxp 



r dr 

n -1 



(5) 



where R = r a A^ 3 . For 6 He, we use the parameter 
set a=0.65 fm, r =1.25 fm, V = -47.4 McV, and 
/ SO =0.93, that reproduces the measured low-energy ri- 
ot phase shifts [3]. We employ the same parameters for 
the density-dependent interaction as in the line 5 of Ta- 
ble II in Rcf. [3]: a nn = -15 fm, E cut = 40 MeV, 
v p = —vq, i? p =2.436 fm, and a p =0.67 fm. The contin- 
uum single-particle spectrum is discretized with a radial 
box of i?box = 30 fm. 

For n Li, we use a=0.67 fm, r =1.27 fm, / so =1.006, 
and Rhox = 40 fm. For the Woods-Saxon potential, a 
deep potential Vb = —47.5 MeV is used for the even par- 
ity states, while a shallow potential Vb = —35.366 MeV 
is adopted for the odd parity states in order to increase 
the s-wave component of the ground state wave function 
[3]. Similar potentials are used also in Ref. [4]. For the 
density dependent force, we use a similar parameter set 
dnn = —15 fm, E cut = 30 MeV, v p = — v 0} R p = 2.935 fm, 
and a p =0.67 fm to that of the line 5 in Table IV in Ref. 
[3] except the value of the energy cutoff. 

For 24 0, we use V = -43.2 MeV, a=0.67 fm, r =1.25 
fm, / SO =0.73, and i?box = 30 fm so that the bound S1/2 
and d 5 / 2 states have empirical single-particle energies 
-2.739 and -3.806 MeV observed in 23 and 21 0, re- 
spectively. For the pairing interaction, we use a nn = —15 
fm, E cut = 30 MeV, v p = 814.2 MeV fm 3 , R p = R, and 
a p =0.67 fm which reproduce the two neutron separation 
energy of 24 0, S 2 „=6.452 MeV. 

The calculated ground state properties are summarized 
in Table 1, where 

(rL) = (* 9 s\(n-r 2 ) 2 \* gs ), (6) 

is the mean square distance between the valence neu- 
trons, and 

<^JH* 3S |(ri+r 2 ) 2 /4|* 3S >, (7) 

is the mean square distance of their center of mass with 
respect to the core. 



TABLE I: Ground state properties of 6 He, n Li, and 24 ob- 
tained with the three-body model with the density-dependent 
delta interaction. The result for 6 He is the same as the line 5 
in Table II in Ref. [3]. 



nucleus 


S2n (r£„) (r^L 2 n) 


dominant 


fraction 




(MeV) (fm 2 ) (fm 2 ) 


configuration 


(%) 


6 He 


0.975 21.3 13.2 


(PS/2) 2 


83.0 


"Li 


0.295 41.4 26.3 


(Pl/2) 2 


59.1 


24q 


6.452 35.2 10.97 


(Sl/2) 2 


93.6 



III. DISCUSSIONS 
A. Ground state properties 

Let us now discuss the spatial correlation of the valence 
neutrons in the ground state, and its influence to the 
dipole excitations near the neutron threshold. To this 
end, we first plot the two-particle density. It is given as 
a function of two radial coordinates, r\ and r 2 , for the 
valence neutrons, and the angle between them, #i 2 . The 
two-particle density can be decomposed into the 5 = 
and 5 = 1 components in the L5-coupling scheme, i.e., 



02(ri,r 2 , 0i 2 ) = P2 (r u r 2 , 12 ) + p 2 (n, r 2 , 6 12 ). (8) 
The explicit expression for the each component is given 

by [1] 



p!=°(ri,r2,0i2) 



EEE 



U'L I I V L 



L 1,3 I ,3 \ / 

x(-) 1 



/2j + l /2j' + l 



21 + 1 V 21' + 1 



(9) 



/of-^ri,^,^) = 



EEE 



U'L I V L\ I V L 



8tt V'H'r^ V^f 0001-10 

L 1,3 I ,3 V / V 

x $o ( r i> r 2) $ i' j' {ri,r 2 )Y LO (0i 2 ) 



2.7 + 1 j n 2j' + 1 
21' + 1' 



(10) 



where I = \/2l + 1. Here, &ij(r, r') is the radial part of 
the two-particle wave function defined as 



n ^ iV /2(l + W) 
x {4>nij{r)(j) n aj{r') + <f> n ij{r')<t>n'ij(r)),(U) 

where n and n' are the radial quantum numbers and 
ot-nn'lj is the expansion coefficient. (f> n ij(r) is the radial 
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FIG. 1: (Color online) The two-particle density for 6 He as 
a function of r\ = ra = r and the angle between the valence 
neutrons, #12. The two-particle density is weighted with a fac- 
tor 4nr 2 ■ 2nr 2 sin #12. The top panel shows the total density, 
while the middle and the bottom panels show the S=0 and 
the S=l components in the LS-coupling scheme, respectively. 



part of the Woods-Saxon single particle wave function. 
Notice that the two-particle density is normalized as 



poo poo /*n 

/ ATrrjdri / r\dr 2 / 27rsin0i 2 dOu P2(»T> ?"2, 0i2, 
Jo Jo Jo 

(12) 

Figures 1,2, and 3 show the (total) two-particle den- 
sity (the top panels) for the 6 He, 11 Li, and 24 nuclei, 
respectively, and their spin decompositions (the middle 
and the bottom panels). These are plotted as a func- 
tion of the radius r\ = r 2 = r and the angle #12, and 
with a weight of Airr 2 ■ 2irr 2 sin #12. As has been pointed 
out in Refs. [1, 6, 12], one observes two peaks in the 
two-particle densities, although the two peaked structure 
is somewhat smeared in 24 O. The peaks at smaller and 
larger #12 are referred to as "di-neutron" and "cigar-like" 
configurations in Refs. [6, 12], respectively. We see that 
the di-neutron part of the two-particle density has a long 
radial tail in 6 He and 11 Li, and thus can be interpreted 
as a halo structure. In contrast, the cigar-like configu- 
ration has a rather compact radial shape. For 24 0, the 
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FIG. 2: (Color online) Same as fig. 1, but for 11 Li. 



di-neutron and the cigar-like configurations behave simi- 
larly as a function of r, and do not show a halo structure. 
Evidently, a large rms radius of 24 is attributed to the 
dominant s-wave component in the ground state wave 
function, rather than the halo effect (see below). 

We find that the spin structure of the two-particle den- 
sity is considerably different among the three nuclei stud- 
ied. In order to see this transparently, we introduce the 
angular density p(0i 2 ) by integrating the radial coordi- 
nates in the two-particle density, i.e., 



p(6i 2 ) = 4tt 



p CO nC 

/ r\dr x 
Jo Jo 



rldr2P2(ri,r 2 ,0i 2 ). (13) 



The angular density is normalized to unity as 



2tt / sin(9i2d6)i2p((9i2) = 1. 



(14) 



Fig. 4 shows the angular density for the 6 He, n Li, and 
24 (with a weight of 27rsin(9i2). The solid line is the 
total density, while the dashed and the dotted lines are 
for the S — and the S — 1 components, respectively. 
The expectation value of the angle 6*12 is 66.33, 65.29, and 
82.37 degree, for 6 He, n Li, and 24 0, respectively. For the 
borromean nuclei 6 He and n Li, the S=0 wave function 
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FIG. 4: The angular density defined by Eq. (13) (weighted 
with a factor 27rsin#) for 6 He, n Li, and 24 0. The solid line 
is for the total density, while the dashed and the dotted lines 
are for the 5=0 and the 5=1 components in the iS-coupling 
scheme, respectively. 



FIG. 3: (Color online) Same as fig. 1, but for 24 0. 



dominates the di-neutron part of the two-particle density. 
In contrast, the cigar-like part has a large 5=1 part in 
n Li, but still the 5=0 component dominates in 6 He. For 
the 24 nucleus, there is no clear separation between the 
di-neutron and the cigar-like type structures. The wave 
function shows a strong correlation typical in the BCS 
type wave function. In fact, the calculated rms radius 
for 24 0, 4.45 fm, is close to that of the bound 2si/ 2 state, 
4.65 fm. A small difference in the rms radii is due to the 
anti-halo effect, caused by the pairing correlation [13]. 

The main features of the angular dependence of the 
two-particle density shown in Fig. 4 can be understood 
in the following way. From Eqs. (9) and (10), one can 
obtain by inserting the values of 3j-symbols that 



P' 17 ! 1 

P S=1 



2^/5 

15 

5 



s ""*--) cx -y o(M + ^^20(M °c cos 2 0^15) 
12) oc -loo (#12) 



as 



^20(^12) cx sin 2 0^16) 



for the configurations — =P3/ 2 or — 

=Pi/2 • When weighted by sin#i 2 , Eq. (15) has 
a peak at #i 2 =35.26 degree and 144.74 degree, while Eq. 
(16) has the maximum at #12=90 degree. This is indeed 
the case for the 6 He nucleus. For the 11 Li, the admixture 
of the (si/2) 2 and (d5/2) 2 configurations perturb this pic- 



ture, and a peak at $12=144.74 degree in the 5=0 com- 
ponent disappears to a large extent. For the 24 nucleus, 
the 5=1 component is largely suppressed since the pure 
(si/2) 2 state cannot form the 5=1 configuration. Also, 
the two-particle density for the pure (si/ 2 ) 2 configuration 
is proportional to |loo| 2 , and thus has a peak at #12 =90 
degree when it is weighted by sin 6*12 . 



B. Dipole excitations 

We next discuss the response of the ground state to 
the dipole field, 

Dm = -je (riF 1M (ri) + r 2 Y 1M (r 2 )). (17) 

Since we obtain the excited 1 _ states by the matrix diag- 
onalization, they appear as discrete states. We smear the 
discrete strength distribution with a smearing function as 



B(E1) = 



r 



tt (E~E k y + r< 



■B k (El), 



(18) 



where Bk(El) is the B(E1) strength for the fc-th excited 
state, 



S Jb (£?l)=3|(*J_|Uo|* i 



(19) 
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FIG. 5: The B(E1) distribution for 6 He, n Li, and 24 0. The 
solid and the dashed lines are obtained with a smearing pro- 
cedure with T = 0.2 and 0.5 MeV, respectively. 



Fig. 5 shows the B(E1) distributions for the 6 He, 11 Li, 
and 24 O. The solid and the dashed lines are obtained 
with the smearing function (18) with T=0.2 and 0.5 MeV, 
respectively. The discrete distributions are also shown. 
The total B{El) strength, J2k B k( E1 ), is 1.31, 1.76 and 
0.97 e 2 fm 2 for the 6 He, n Li, and 24 0, respectively. 

We notice that the strong threshold peak appears in 
the response of 6 He and n Li. In the cluster model within 
the plane wave approximation, the peak of the strength 
function appears at 1.65c, where 5 C is the cluster sep- 
aration energy [14-16]. The calculated peaks in Fig. 5 
are at 1.55 and 0.66 MeV for 6 He and 11 Li, respectively. 
These peaks are very close to 1.6 times the two neutron 
separation energy, 1.66 , 2 r( =l-56 MeV for 6 He and 0.47 
MeV for n Li. This similarity suggests the existence of 
strong di-ncutron correlations in these nuclei. A small 
difference between the peak energy in Fig. 5 and 1.652« 
for n Li is due to the large configuration mixing of S1/2 
state (22.7%) in the ground state. In contrast, for 24 0, 
the peak (4.78 MeV) is below the two neutron separation 
energy, and is rather close to 1.6 times single particle en- 
ergy for the 2s 1/2 state, that is, 1.6x 2.74 = 4.38 MeV. 
Therefore, the di-neutron correlation does not seem to 
play a major role in the dipole response in this nucleus. 

In order to see the threshold effect more clearly, we plot 
in Fig. 6 the transition density for the strong El peaks 
for each nucleus. As a comparison, we also show the 
5=0 component of the transition density by the dashed 
line. For 6 He and 11 Li, the transition density shows a 




r (fm) 

FIG. 6: The transition density for the state close to the peak 
in the the B(E1) distribution. The top, middle, and bottom 
panels are for 6 He, n Li, and 24 0, respectively. The solid 
line indicates the total density, while the S=0 component is 
denoted by the dashed line. 



nodal structure and changes its sign, that is typical in 
the coupling to the continuum spectrum. Also, the 5=0 
component plays a significant role, supporting the impor- 
tance of the di-neutron correlation in these nuclei. On the 
other hand, such clear nodal structure is not seen in the 
transition density of 24 0. The transition density seems 
to consist of a coherent sum of the 5=0 and the 5=1 
components. The dipole excitation in 24 is therefore 
interpreted as a coherent superposition of particle-hole 
excitations, rather than the continuum excitations, as in 
stable nuclei where the continuum effect is much less im- 
portant. 

We mention that, while the Coulomb breakup for the 
one-neutron halo 11 Be is now well established, that for 
the two-neutron halo 11 Li is still in dispute showing large 
discrepancies between the experimental data taken by 
three different groups [17]. Recently, Nakamura and his 
collaborators performed the Coulomb dissociation exper- 
iments of n Li on 208 Pb target with much higher statis- 
tics and with much less ambiguities caused by cross talk 
events in detecting two neutrons [18]. They observed a 
sharp peak at E oxp ~0.6MeV and the integrated strength 
B cxp (El)=1.5±0.1e 2 fm 2 for E< 3.3MeV, which are con- 
sistent with the present results E pea k=0.66 MeV with 
the calculated strength B cal (El)=1.31 e 2 fm 2 for E< 
3.3MeV. Similar calculated results were also reported in 
Rcf. [3]. 
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IV. SUMMARY 

We have studied the role of di-neutron correlations in 
weakly bound nuclei on the neutron drip line using a 
three-body model. The model Hamiltonian consists of a 
Woods-Saxon potential between a valence neutron and 
the core, and of a density-dependent pairing interaction 
among the valence neutrons. We applied this model to 
the borromean nuclei, 6 He and 11 Li, as well as a skin 
nucleus 24 0. For the borromean nuclei, 6 He and 11 Li, 
we found that the two-particle density has a two peaked 
structure, one peak at a small opening angle between 
the valence neutrons and the other at a large angle (the 
'di-neutron' and 'cigar- like' configurations, respectively). 
We found that the former is dominated by the S—0 con- 
figurations in the LS coupling scheme and has a long 
tail. On the other hand, the latter has a compact shape, 
and is dominated by the S=l configuration for 11 Li while 
by the 5=0 configuration for 6 He. For 24 0, there is no 
clear separation between the di-neutron and the cigar- 
like configurations, and the ground state is dominated 
by the S=0 configuration. 

We have also studied the dipole response of these nu- 
clei within the same model. We found strong threshold 



peaks in the response of 6 He and n Li nuclei, where the 
transition density shows the importance of the coupling 
to the continuum. The S=0 configuration, thus the di- 
neutron correlation, was found to have a large contribu- 
tion to the transition density for these peaks. On the 
other hand, no clear sign of the continuum coupling was 
seen in the response of 24 0. The transition density for the 
low energy dipole strength of 24 O consists of a coherent 
sum of the S—0 and the S=l components, and therefore 
the di-neutron correlation plays a much less important 
role in 24 than in the borromean nuclei. 

Recently, a new Coulomb breakup measurement for 
11 Li has been undertaken at RIKEN [18]. It would be 
interesting to perform a similar experiment also for the 
24 O nucleus and see a difference between them as we 
discussed in this paper. 
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